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ABSTRACT 

We examine the test-particle solution for diffusive shock acceleration, based 
on simple models for thermal leakage injection and Alfvenic drift. The critical 
injection rate, £ c , above which the cosmic ray (CR) pressure becomes dynamically 
significant, depends mainly on the sonic shock Mach number, M, and preshock 
gas temperature, 7\. In the hot-phase interstellar medium (ISM) and intracluster 
medium, £ c < 1(T 3 for shocks with M < 5, while £ c « l(T 4 (Ti/10 6 K) 1/2 for 
shocks with M > 10. For T\ = 10 6 K, for example, the test-particle solution 
would be valid if the injection momentum, p in j > 3.8p t h (where p t h is thermal 
momentum). This leads to the postshock CR pressure less than 10% of the 
shock ram pressure. If the Alfven speed is comparable to the sound speed in 
the preshock flow, as in the hot-phase ISM, the power-law slope of CR spectrum 
can be significantly softer than the canonical test-particle slope. Then the CR 
spectrum at the shock can be approximated by the revised test-particle power- 
law with an exponential cutoff at the highest accelerated momentum, p max (t). 
An analytic form of the exponential cutoff is also suggested. 

Subject headings: acceleration of particles — cosmic rays — shock waves 



1. Introduction 

Suprathermal particles are produced as an inevitable consequence of the formation 
of collisionless shocks in tenuous astrophysical plasmas and they can be further acceler- 
ated to very high energies through int eractions with resonantly scattering Alfven waves 



in the converging flow across a shock (jBelll Il978l ; iDruryl Il983l ; iBlandford &: Eichlerl 11987 



Malkov fc Druryll200ll ). The most attractive feature of the diffusive shock acceleration (DSA) 
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theory is the simple prediction of the power-law momentum distribution of cosmic rays (CRs), 
f(p) oc p _3ct /( <t_1 ) (where a is the shock compression ratio) in the test particle regime. For 
strong, adiabatic gas shocks, this gives a power-law index of 4, which is reasonably close to 
the observed, 'universal' index of the CR spectra in many environments. 

The nonthermal particle injection and ensuing acceleration at shocks depend mainly 
upon the shock Mach number, field obliquity angle, and the strength of the Alfven turbu- 
lence responsible for scattering. At quasi-parallel shocks, the shock Mach number is the 
primary parameter that determines the CR acceleration efficiency, while the injection frac- 
tion, £ (the ratio of CR particles to the total particles passed through the shock), is the 
secondary parameter. Detailed nonlinear treatments of DSA predict that at strong shocks, 
with a small fraction of £ > 10 -4 , a significant fraction of the shock kinetic energy is trans- 
ferre d to CRs and there are highly nonlinear ba ck-reactions from CRs to the underlying 



flow (IBerezhko fc Volk 1119971 ; iKang fc Jonesll2007l ). Indeed, multi-band observations of non 



thermal radio to 7-ray emissions from several supernova remnants (SNRs) have been sue- 
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It has been recognized, however, that the CR spectrum at sources, N(E), predicted for 
shocks strongly modified by CR feedback may be too flat t o be c onsistent with the observed 
flux of CR nuclei at Earth, J(E). Recently lAve et al.l (120091 ) analyzed the spectrum of 
CR nuclei up to ~ 10 14 eV measured by TRACER instruments and found that the CR 
spectra at Earth can be fitted by a single power law of J(E) oc E~ 2 ' 67 . Assuming an energy- 
dependent propagation path length (A oc E~~ - 6 ), they suggested that a soft source spectrum, 
N(E) oc E~ s with s ~ 2.3 — 2.4, is preferred by the observed data. This is much softer than 
the CR spectrum that the nonlinear DSA predicts for strong SNRs, which are believed to be 
the main accelerators for Galactic CRs up to the knee energy around 10 15,5 eV. Thus, in order 
to reconcile the DSA prediction with the observed J(E), the bulk of Galactic CRs should 
originate from SNRs in which the CR acceleration efficiency is 10 % or so (i.e., roughly in 
the test-particle regime). Such inefficient acceleration could be possible for SNRs in the hot 
phase of the interstellar medium (ISM) (i.e., low shock Mach number shocks) and for the 
inject fraction smaller than 10 -4 (jKangil2010l ). 



The scattering by Alfven waves tends to isotropize the CR distribution in t he wave 
fram e, which may drift upstream at Alfven speed with respect to the bulk plasma (ISkilling 
19751 ). This Alfvenic drift in the upstream region reduces the velocity jump that the par- 
ticles experience across the shock, which in turn softe ns the CR spectrum beyond th e 
canonical test-particle slope (s = 2 for strong shocks) (IKang! |2010| ; ICaprioli et al.ll2010l ). 
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Moreover, the Alfvenic drift in amplified magnetic fields both upstream and downstream 
can drastically soften the accelerated particle spec trum even in nonlinear modified shocks 
rtzirakashvili and Ptuskm1l2008l : IPtuskin et o/lboioh . 



At collisionless shocks suprathermal particles moving faster than the postshock thermal 
distribution may swim through th e MHD waves and leak upstream across the shocks a nd get 
injected into the CR population rtMalkovl [l~998l ; bieseler et allboool : kang et al.ll2002l ). But 
it is not yet possible to make precise quantitative predictions for the injection process from 
first principles, because complex plasma in teractions among CRs, waves, and the underlying 



gas flow are not fully understood yet (e.g., iMalkov fc Druryll200ll ). Until plasma simulations 



such as hybrid or particle-in-cell simulations reach the stage where the full problem can be 
treated with practical computational resources, in the studies of DSA we have to adopt a 
phenomenological injection scheme that can emulate the injection process. 

In this paper, we wi l l exam ine the relation between the thermal leakage injection model 
described in iKang et al.l ( 120021 ) and the time-dependent test-particle solutions for DSA. The 
basic models are described in §2, while the analytic expression for the CR spectrum in the 
test-particle limit is suggested in §3. Finally, a brief summary will be given in §4. 



2. BASIC MODELS 



In the kinetic DSA approach, the following diffusion-convection equation for the pitch- 
angle- aver aged distribution function, f(x,p,t), is solved along with suitably modified gasdy- 
namic equations: 



dt 



— - [U + u w)t^ = - 



p d{u + u h 
3 dx 



dp dx 



n(x,p)- 



where n(x,p) is the spatial diffusion coefficient and u m is the drift speed of the local Alfvenic 
wave turbulence with respect to the plasma (lSkillinglll975f ). We consider only the proton 
CR component. 



2.1. Alfvenic Drift Effect 

Since the Alfven waves upstream of the subshock are expected to be established by the 
streaming instability, the wave speed is set there to be u w = —v^. Downstream, it is likely 
that the Alfvenic turbulence is nearly isotropic, hence u w = there. As a result, the velocity 
jump across the shock is reduced, and the slope of test-particle plower-law spectrum should 
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be revised as 

3(m-^) 3a(l - M/) 

where «i and u<i are the upstream and downstream speed, respectively, in the shock rest 
frame, o = U\jui = pij p\ is the shock compression ratio, and va and Ma = ui/va are the 
Alfven speed upstream and Alfven Mach number. Hereafter, we use the subscripts '1', and 
'2' to denote conditions upstream and downstream of the shock, respectively. Thus the CR 
spectrum would be softer than the canonical power- law spectrum with the slope, 3er/(er — 1), 
unless M A > 1. 

The left panel of Figure 1 shows the revised test-particle slope function of the 

sonic Mach number, M, for different Alfven speeds, v a = 5 • c s (where c s is the upstream 
sound speed). In the hot-phase ISM of T rj 10 6 K with the hydrogen number density n# ~ 

0. 003 cm -3 and the magnetic field strength B w 5/xG, the sound speed is c s ~ 150 km s _1 and 
the Alfven speed is va ~ 170 km s -1 . So 5 ~ 1 is a representative value. If 8 ~ {Pb/ P g ) 1 ^ 2 ~ 

1, the Alfven drift effect is significant for Alfven Mach number, Ma ~ M < 30. Consequently, 
this effect reduces the CR acceleration efficiency. Of course, it is not important for strong 
shocks with u s ^> c s ~ t>A {i.e., Ma > 30). 



2.2. Thermal Leakage Injection Model 



Since the velocity distribution of suprathermal particles is not isotropic in the shock 
frame, the diffusion-convection equation cannot directly follow the injection from the non- 
diffusive thermal pool into the diffusive CR popul ation. Here we adopt the thermal leakage 
injection model that was originally formulated by iGieseler et al.l ( 120001 ) based on the calcu- 
lations of iMalkovl ( 119981 ). In this model particles above a certain injection momentum p in ^ 
cross the shock and get injected to the CR population. We adopt a smooth "transparency 
function", r esc (e j B,f), that expresses the probability of suprathermal particles at a given ve- 
locity, v, leaking upstream through the postshock MHD waves. One free parameter controls 
this function; €b = E>o/B±, the ratio of the general magnetic field along the shock normal, 
Bq, to the amplitude of the postshock MHD wave turbulence, B i . Although plasma h ybrid 
simulations and theories both suggested that 0.25 < e# < 0.35 OMalkov fc Volklll998l ). the 
physical range of this parameter remains to be rather uncertain due to lack of full under- 
standing of relevant plasma interactions. Since r esc increases gradually from zero to one in 
the thermal tail distribution, the "effective" injection momentum can be approximated by 



1.07, 



p inj w 1.17m p w 2 (l + - — ) = Qmj{M,e B )p, 



th 



(3) 
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where pth = \f^rn^B^2 is the thermal peak momentu m of the immedia te postshock gas 
with temperature T2 and ks is the Boltzmann constant ( Kang et al.lhooj ). 



The right panel of Figure 1 shows the value of Qinj as a function of M for three values of 
€b = 0.21, 0.23, and 0.27, which represents "inefficient", "moderately efficient", "efficient" 
injection cases, respectively (see Fig. 4 below). At weaker shocks the compression is smaller 
and so the ratio Uiju\ is larger. For stronger turbulence (larger B±, smaller es) it is harder 
for particles to swim across the shock. So for both of these has to be larger. Hence 

the value of Qi n j(M, eb) is larger for weaker shocks and for smaller eg, which leads to a lower 
injection fraction. 

In our thermal leakage injection model, the CR distribution function at p in j is then 
anchored to the postshock Maxwellian distribution as, 

/inj = /(Pmj) = ^ Pth exp(-gf nj ), (4) 

where n 2 is the postshock proton number density and the distribution function is defined 
in general as J Airp 2 j (p)dp = n. For the test-particle power-law spectrum, the value of 
Qmj determines the amplitude of the subsequent suprathermal power-law distribution as 
f(p) = /inj ■ (p / Pmj)~ qtp ■ Then the CR injection fraction can be defined as 

which depends only on the ratio Qi n j and the slope q tp , but not on the postshock temperature 
T 2 . For Qinj = 3.8, for example, £ = 6.6 x 10~ 5 /(g tp — 3), which becomes £ = 6.6 x 10~ 5 for 
strong shocks with q tp = 4.0. 



2.3. Bohm-type Diffusion Model 



In modeling DSA, it is commonly assumed that the particles are resonantly scattered 
by self-generated waves, so the Bohm diffusion model can represent a saturated wave spec- 



trum (i.e., the mean scattering length, A 



91 



where r g is the gyro-radius). Here, we 



adopt a Bohm-type diffusion coefficient that includes a weaker non-relativistic momentum 
dependence, 



k(x,p) 



K 



P 



m p c 



p(x) 
L Pi 



(6) 



where the coefficient k* = m p c 3 / (3eBo) depends on the upstream mean field strength. The 
case with m = 1 approximately accounts for the compressive amplification of Alfven waves. 
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The mean acceleration time for a particle to reach p max from pi n j in the test-particle 
limit of DSA theory can be approximated by 



u%—v A - u 2 



Pmax 



Ki K 2 \ dp 

- va u 2 J p ' 



(7) 



if we assume the bulk drift of waves with v A in the upstream region (e.g., iDruryl LL983J ). 
Then the maximum momentum can be estimated by setting t = t acc as 



3a[l + (1 - M/)ct 



X-m] 



K 



K 



where u s = U\ is the shock speed (IKang et al.l 120091 ). For the case of m = 1, the typical 
value of the parameter, f c = a(l — M^)(a — 1 — aM^ 1 )/{3a[l + (1 — M A 1 )a 1_m ]}, is ~ 1/8 
in the limit of M A > 1 and M > 1. 



3. TEST-PARTICLE SPECTRUM 

If the injection is inefficient, especially at weak shocks, the CR pressure remains dy- 
namically insignificant and the test-particle solution is valid. Caprioli et al. (2009) (CBA09 
hereafter) derived the analytic solution for a steady- state, test-particle shock with a free- 
escape boundary (FEB) at a distance £feb upstream of the shock (i.e., f(x > xfeb) = 0). 
For a diffusion coefficient that depends on the momentum as n(p) = K*(p/m p c) a , the CR 
distribution at the shock location, x s , is given by 

ft P (x s ,p) = fo ■ exp 

where z = p/p*, z in j = Pmj/p*, fo = fmj, and p*/m p c = (xfebw s / ' k*) 1 / 01 is the cutoff momen- 
tum set by the FEB. This expression can be re-written as, 

ft P (x s ,p) = / inj • • exp [-q tp C(z)\ , (10) 

Pinj 

where the function C(z) is given by 

r dz' i 

We show the function C(z) for a = 0.5 and 1 in the left panel of Figure 2. For 2 C 1, C(z) 
is small and so exp [—q tp C(z)} = 1, as expected. For z 1, C(z) m z a = (p/p*) a . But this 



-Qt P 



dz' 



z' 1 - exp(-l/z /Q ) 



(9) 
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regime (p ^> p*) is not really relevant, because the resulting ft p (x s ,p) is extremely small. 
We are more interested in the exponential cutoff where p ~ p*. Figure 2 shows that C(z) 
increases much faster than z a near z ~ 1. In fact, at z ~ 1, approximately C(;2) ~ 0.29z 2 
for a = 1 and C(z) m 0.58z for a = 1/2. Thus equation ( flOi) can be approximated by 

0.29g tp p 2q " 



ftp(x a ,p) ~/inj- (— ) gtp -exp 

Pinj 



a p* 



(12) 



Kang et al.l (120091 ) showed that the shock structure and the CR spectrum of time- 
dependent, CR modified shocks with ever increasing p max (t) are similar to those of steady- 
state shocks with particles escaping through the upper momentum boundary, i.e., f(p > 
Pub) = 0, if compared when p max (£) = Pub (see their Figs. 10-11). They also showed that the 
exponential cutoff in the form of exp[— k(p/p mauX ) 2a ] matches well the DSA simulation results 
for CR modified shocks. In the same spirit, we suggest that equation ffTO]) could represent the 
CR spectrum at the shock location for time-dependent, test-particle shocks without particle 
escape, in which the cutoff momentum is determined by the shock age as in equation (jHJ), 
i.e., p* ~ p max (t). 

The distribution function f(x, p max ) in the upstream region decreases roughly as exp[— x/ld(p m &x)}, 
where the diffusion length for p max is 

^(Pmax) = ^ PmaX ^ = f C U s t. (13) 

CBA09 spectrum in equation (TTU]) was derived from the FEB condition of f(x > £feb,p) = 
for steady-state shocks, while f(x,p) — > only at x — > oo (upstream infinity) for time- 
evolving shocks without particle escape. So we presume that the cutoff momentum can be 
found by setting the location of FEB at xfeb = C ' Uftnax), where ( ~ 1. From the condition 
that p*/m p c = ((ld(Pm ax )u s /K*) 1/a , we find p* = ( ■ p max . 

The right panel of Figure 2 shows the test-particle solution from a time-dependent DSA 
simulation, in which the dynamical feedback of the CR pressure was turned off. Contrary 
to CBA09 case, no FEB is enforced in this simulation, so the shock does not approach to a 
steady state, but instead evolves in time. As the CRs are accelerated to ever high energies 
(Pmax oc t), the scale length of the CR pressure increases linearly with time, /d(Pmax) oc u s t. 
So the shock structu re evolves in a self-s imilar fashion, depending only on the similarity 



variable, x/(u s t) (see iKang fc Jonesl 120071 ). By setting p* = 1.2p max (t) (i.e., ( = 1.2) and 
also by adopting the value of / in j from the DSA simulation result, we calculated ft P (x s ,p) 
according to equation (TTU]) . As can be seen in the figure, the agreement between the numerical 
DSA results and the analytic approximation is excellent. Thus we take equation ([101) as the 
test-particle spectrum from DSA, where q tp , Pmj, fmj, and p* ~ 1.2p max (t) are given by 
equations (T2]), ([3]), (TjJ, and ([8]), respectively. 
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Figure 3 shows some examples of the test-particle spectrum given in equation ffTO]) . We 
consider the shocks propagating into the hot-phase of the ISM of T\ = 10 6 K or a typical 
intracluster medium (ICM) of T\ = 10 7 K. The shock speed is given by u s = M ■ c s , where 
the sound speed is c s = 150 km s _1 (T 1 /10 6 K) 1 / 2 . For all the cases, we assume a constant 
cutoff momentum, p* = 10 6 GeV/c, which is close to the knee energy in the Galactic cosmic 
ray spectrum. For typical hot-phase ISM, 5 = va/c s ~ 1 as mentioned before. For typical 
ICM, tih ~ 10~ 3 cm~ 3 and B m 1 — 5//G, so 5 ~ 0.5 is taken here. For typical test-particle 
limit solutions, we adopt e# = 0.21 to specify p^ given in equation fl3]), which determines 
the anchoring point where the test-particle power-law begins. This choice of results in 
the injection rate £ < 10~ 4 and the postshock CR pressure P C)2 /(piwf) < 0.1. As can be 
seen in Figure 3, for stronger (faster) shocks, the postshock gas is hotter, the amplitude /; n j 
is higher and the power-law spectrum is harder. 

Then the CR pressure at the shock position can be calculated by 

, s 47r f°° , / \ P 4 dp , . 

For strong shocks with q tp = 4, with the test-particle spectrum in equation ( !T0|) . P c oc 
/injPfnj ki(p*/m p c). Then, with a constant cutoff p*, P c oc exp(— Qf^Qf^Pth- So for a fixed 
value of Qi Q j (or fixed injection fraction £), P c oc p t h oc u s . Figure 4 shows the fraction of 
injected particles and the postshock CR pressure calculated by adopting the test-particle 
spectrum given in equation (fTOl . The same p* = 10 6 GeV/c is chosen as in Figure 3. The 
quantities, ra cr 2 and P C)2 do not depend sensitively on the assumed value of p* for weak 
shocks, since the power-slope qt p is greater than 4. But for strong shocks (M > 30) where 
q tp ~ 4 (see Fig. 1), the CR pressure increases logarithmically as P c oc \n(p* / m p c) . Several 
values of T 1; t B (or Qinj), and 5 = va/c s are considered. In general, for fixed values of ee (or 
Qinj) and 5, the ratio P c ,2/ '{pi u2 s ) increase strongly with the shock Mach number for shocks 
with M < 10, because of the strong dependence of £ (or Qmj) on M for weaker shocks. 
But for shocks with M > 10, £ becomes independent of M and so P c discussed 
above. So the CR pressure relative to the shock ram pressure, P c ^/{piU 2 s ) oc u^ 1 , that 
is, it becomes smaller at faster shocks. Of course, in the nonlinear DSA regi me, the ratio 



Es^j/ jpiUg) increases with the shock Mach number and saturates at about 1/2 flKang et al. 



20091 ). 



The top panels of Figure 4 show how the CR pressure depends on and 5. For a given 
Mach number, the CR pressure increases strongly with €b, because of the exp(— Qf n $) factor. 
Obviously, the CR pressure becomes smaller for larger 5 because of softer power-law spectra 
at weaker shocks with M < 30. For e B = 0.21 and 6 = 1, £ < 10~ 4 and P C)2 /(pi^) < 0.1, so 
the test-particle solution would provide a good approximation. For e# = 0.23, on the other 
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hand, the injection fraction becomes £ ~ 10 4 — 10 3 , and the test-particle solution is no 
longer valid for M > 5. For weak cosmological s hocks with M < 3, typically found in the 



hot ICM (e.g., iRyu et a/.ll2003l ; iKang et aZ.ll200?t ) . even for a rather large value of e# = 0.27 



the injection fraction is smaller than 10~ 3 and P c ,il 'piu 2 < 0.01 So we can safely adopt the 
test-particle solution for those weak shocks, unless there are abundant pre-existing CRs in 
the preshock flow. 

The middle panels show the cases with the same Qmj, independent of M. For these 
cases, T\ = 10 6 K, 5 = 1, and p max = 10 6 GeV/c. With the same Q in j, the injection fraction is 
almost independent of M except for weak shocks with M < 5. For Q in j = 3.8, P c ^/ ' P\u 2 < 0.1 
for all shocks. One can see that Qi n j ~ 3.8 is the critical value, above which the injection 
fraction becomes £ < 10~ 4 and the ratio P c ^/(piu 2 ) < 0.1. Hence, if p in j > 3.8p t h, the CR 
injection fraction is small enough to guarantee the validity of test-particle solution. But once 
again one should note that P c oc Inp* for strong shocks. 

The bottom panels show the cases in which the preshock temperature is T\ = 10 5 — 10 7 K. 
Since the ratio P c ^/(piU 2 ) oc ^u^ 1 and £ does not depend on T\, P c ^/(p\U 2 ) oc £T 1 _1 ' /2 for a 
given Mach number, M = u s /c s . So we chose e# ~ 0.20 — 0.22 for different T 1; which results 
in £ ~ 10 _4 (Ti/10 6 K) 1//2 . This gives the similar value of P c ,2/ (p\u 2 s ) ~ 0.1 for three values 
of T\. For these shocks, the test-particle solution would be valid. 

When P c 2 /(piU 2 ) > 0.1, the nonlinear feedback of the diffusive CR pressure becomes 
important and the evolution of CR modified shocks should be followed by DSA simulations. 
Figure 5 compares the evolution of a slightly modified M = 5 shock (eb = 0.27) with 
that of a test-particle shock (e# = 0.2). In the CR modified shock, the upstream flow 
is decelerated in the precursor before it enters the gas subshock. So the quantities at far 
upstream, immediately upstream and downstream of the subshock are subscripted with '0', 
T', and '2', respectively. For the test-particle shock, p 1 = p and 7\ = T . Here T = 10 6 K 
and Va/c s = 0.42. The simulations start with a purely gasdynamic shock at rest at x = 0, 
initialized according to Rankine-Hugoniot relations with Uq = — 1, po = 1 and a gas adiabatic 
index, 7 5 = 5/3. There are no pre-existing CRs. 

The test-particle spectrum given in equation ffTOl) with p* = 1.2p max at t/t = 10 is also 
shown for comparison (dot-dashed lines) in the bottom panels. In the test-particle shock 
with tB = 0.2, both P Ct 2/(pou 2 ) ~ 0.005 and f(x s ) from the DSA simulation agree well with 
the test-particle solution given in equation (flOl) . as expected. 

If we were to take the test-particle spectrum with tB = 0.27, we would obtain £ = 
1.74 x 10 -3 and P C 2 / (piu 2 ) = 1.17, which is unphysical. In the CR modified solution 
from the DSA simulation, however, £ w 3.6 x 10~ 4 and P c ,2/(piu 2 ) ~ 0.1. The postshock 
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temperature T2 is reduced about 17 % in the CR modified solution (due to higher pi and 
lower p 9t 2), compared to that in the test particle solution. But u 2 and so p in j remain about 
the same. As a result, the amplitude / in j is lower than that of the test-particle spectrum (see 
the bottom right panel of Fig. 5) and so the injection rate is reduced in the CR modified 
solution. The distribution function f(x s ,p) from the DSA simulation is slightly steeper for 
p/rripC < 10 and slightly flatter for p/m p c > 10 than the test-particle power-law, because 
the flow velocity is slightly modified. This demonstrates that the DSA saturates in the limit 
of efficient injection through the modification of the shock structure (i.e., a precursor plus 
a weak gas subshock), which in turn reduces the injec tion rate. Thus the ratio P c ^/(pxu 2 ^) 



approaches to ~ 1/2 for strongly modified CR shocks ( IKang fc Jonesl 120071 ) 



Finally, we find that the volume integrated spectrum contained in the simulation volume 
can be obtained simply from F(p) = J f(x,p)dx ~ ftp(x s ,p)u 2 t. This provides the total CR 
spectrum accelerated by the age t. 



4. SUMMARY 

Although the nonlinear diffusive shock acceleration (DSA) involves rather complex 
plasma and MHD processes, the test-particle solution may unveil some simple yet essen- 
tial pictures of the theory. In this study, we suggest an analytic form for the CR spectrum 
from DSA in the test-particle regime, based on simple models for thermal leakage injection 
and Alfvenic drift of self-generated resonant waves. 

If the particle diffusion is specified (e.g., Bohm diffusion), the shock Mach number is 
the primary parameter that determines the efficiency of diffusive shock acceleration. For a 
given shock Mach number, the fraction of injected CR particles becomes the next key factor. 
Since the postshock thermal velocity distribution at the injection momentum determines the 
amplitude of the power-law distribution in the thermal leakage injection model, the ratio 
Qinj = Pmj/Pth is the key parameter that controls the CR injection fraction and in turn 
determines the CR acceleration efficiency On the other hand, as a result of the drift of 
Alfven waves in the precursor, the power- law slope should be revised as in equation (J2]), 
which leads to the CR spectrum much steeper than the canonical test-particle power-law. 
This effect is negligible for shocks with the Alfvenic Mach number, M4 > 30. 

For shocks with the sonic Mach number M > 10, depending on the preshock temperature 
Ti, the injection fraction, £ < £ c w 10 _4 (Ti/10 6 K) 1 / 2 would lead to the downstream CR 
pressure, P c ,2/(pxu1) < 0.1. The exact values depend on other parameters such as va- In 
that case, the CR spectrum at the shock location can be described by the test-particle 
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power-law given in equation f lTUj) . in which the amplitude, /i n j, is fixed by the postshock 
thermal distribution at p in j given in equation (jlj). For supernova remnants in the hot-phase 
of the ISM with Ti = 10 6 K, for example, the CR injection fraction becomes less than 10 -4 , 
if Qinj ^ 3.8 (or e B < 0.21). For weaker shocks with M < 5, the test-particle solution is 
valid even for larger injection fraction, so £ c < 10~ 3 . 

We have shown that the CR spectrum at the shock location in time-dependent, test- 
particle shocks without particle escape could be approximated by the analytic solution given 



i n equ ation (flUI) . which was derived for steady-state, test-particle shocks by ICaprioli et al. 



(120091 ). with the cutoff momentum set as p* ~ 1.2p max (t). If the CR injection is inefficient, 
which should be true for weak shocks with M < 5 found in the intracluster medium, the 
test-particle solution presented in this paper should provide a good approximation. Figure 
4 should provide guidance to assess if a shock with specific properties can be treated with 
the test-particle solution. 

With the injection rate greater than £ c , especially for shocks with M > 5, the spectrum 
deviates from the test-particle form due to the modified flow structure caused by the diffusive 
CR pressure. In fact, the DSA efficiency saturates in strongly modified CR shocks, because 
the postshock temperature ge ts lower and so th e injection rate is reduced. Based on the 



results of the DSA simulations, iKang et al.l ( 120091 ) suggested that CR-modified shocks evolve 



self-similarly once the total pressure is dominated by relativistic particles, and that the CR 
spectrum at the subshock can be approximated by the sum of two power laws with the 
slopes determined by the subshock and total compression ratios with an exponential cutoff 

&t Pmax{t)- 
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Fig. 1. — Left: The test-particle power- law slope, gtp, revised by including the Alfvenic drift 
(Eq. [2]), is shown as a function of sonic Mach number for four different values of 5 = va/c s . 
Right: The ratio Q in j = Pmj/Pth is shown for three values of e# = B /Bj_ = 0.21, 0.23 and 
0.27. 
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Fig. 2. — Left: The function C(z) defined in Eq. ( fTTT) is shown as a function of z = p/p* 
for «;(p) = K*(p/m p c) a . Right: The CR distribution at the shock position, / s -p 4 (in units of 
n\m p c), is shown for a Mach 3 shock. The solid line shows the results of a time-dependent 
DSA simulation without particle escape in the test-particle regime (i.e., no CR feedback 
to the flow). The dotted line represent the test-particle spectrum given in Eq. ([TO]) with 

P 1 • 2pmax • 
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Fig. 3. — Left Test-particle spectra given in Eq. (ITOjl with e# = 0.21, p* = 10 6 GeV/c, 
T\ = 10 6 K, and 5 = va/c s = 1.0. Right: Same as the left panel except T\ = 10 7 K, and 
va/c s = 0.5. Each curve is labeled with the shock Mach number, M, and the shock speed is 
u s = M ■ 150 km s'^Tx/lO 6 ^ 1 / 2 . 
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Fig. 4. — The fraction of CR particles, £ = n cr 2 /^2 (left panels), and the downstream CR 
pressure in units of the shock ram pressure (right panels) are shown for the test-particle 
spectrum given in Eq. ( TTU1) with a fixed p* = 10 6 GeV/c. The shock speed is specified by 
u s = M ■ 150 km s _1 (Ti/10 6 K) 1 / 2 . Upper panels: Shocks with the preshock temperature, 
T\ = 10 6 K. Three values of e# = 0.21, 0.23 (long dashed lines), and 0.27 (dot-dashed) are 
considered with 5 = va/c s = 1.0. For the case with e# = 0.21, three cases with 5 = (dashed 
lines), 0.5 (dotted), and 1.0 (solid) are shown. Middle panels: The same model shocks as 
the upper panels except a constant ratio Qmj = Pmj/pth = 3.6 (dashed lines), 3.8 (solid), 4.0 
(dotted), are shown. The Alfen speed is va/c s = 1.0. Lower panels: Shocks propagating 
into different temperature media, 7\ = 10 5 K (with e# = 0.20, dotted lines), 10 6 K (with 
€b = 0.21, solid lines), and 10 7 K (with e# = 0.22, dashed lines) are shown. The Alfen speed 
is va/c s = 1.0 for all three cases. 
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Fig. 5. — Time-dependent DSA simulation results for a M = 5 shock for e# = 0.20 (low 
injection rate, left panels) and e# = 0.27 (high injection rate, right panels). Here T\ = 10 6 K, 
va/c s = 0.42 and k = 10~ 3 (p / 'm p c) (po / p). At the last time epoch (t/t Q = 10) the cutoff 
momentum becomes p* = 1.2p max m 10 3 GeV/c. The shock structure is shown at t/t a = 0.5 
(long dashed lines), 1 (dotted), 5 (dashed), and 10 (solid) as a function of the similarity 
variable, x/{u s t). The (magenta) dot-dashed lines in the bottom panels represent the test- 
particle spectra given in Eq. ( FlOl at t/t Q = 10. Note that the dot-dashed line (analytic 
solution) and solid line (numerical solution) almost coincide with each other in the case with 
€b = 0.2. Here t D is a normalization constant. 



